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ABSTRACT 



Due to the lack of longitudinal stiffening along 
submarine hulls, they are inherently vulnerable to 
underwater explosions directly off the bow or stern. ' 
Accordion-like deformations of the hull are set into motion 
which could cause dynamic amplification in the transient 
response of attached substructures. In underwater shock 
acceptance tests of internal equipment, this interaction is 
created by exploding a charge in a fore and aft configura- 
tion with the submerged shock test vehicle (SSTV) . With the 
increasing availability of large computers and the rapid 
development of numerical methods, several computer codes 
have been written to predict equipment response to 
underwater shocks. Using the ELSHOK (ELASTIC SHO CK) code, 
this investigation studies the effect of hull/substructure 
interaction on stiffened shell response at resonance 
following an end-on load. The transient response of the ■ 
coupled shell/substructure system from tapered and 
conventional charges of equivalent impulse is examined in 
this study. 
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I. 



INTRODUCTION 



A. BACKGROUND 

Submarines are inherently vulnerable to explosions 
detonated directly off the bow or stern. While they have 
adequate protection against side loads with transverse 
bulkheads and frames, there is little longitudinal 
protection along the hull. As a result of this, they are 
susceptible to accordion-like deformations of the hull upon 
impact from end-on loads. It is possible that the frequency 
of the hull motions could excite internal equipment into 
resonance thereby causing unacceptable damage. This 
investigation examines the coupled elastic response of a 
shell/substructure system when subjected to an end-on 
underwater explosion (UNDEX) as the mass and stiffness of 
the internal substructure is varied. 

The current specifications for conducting underwater 
shock tests on submarine-installed equipment are contained 
in MIL-S-901D [Ref. 1]. This document specifies the 
explosive charge weight and geometry of the test, as well as 
the mounting and orientation of the equipment being tested 
within the submerged shock test vehicle (SSTV) . The SSTV is 
nothing more than a ring-stiffened cylindrical shell with 
circular endplates designed to simulate the hull motions 
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created by an underwater shock. Recognizing the potential 
for unacceptable damage to hull mounted equipment from an 
end-on load, the specification requires the first of four 
UNDEX's to be performed in a fore and aft configuration. 

The suitability of an equipment design or installation is 
evaluated according to its ability to function as intended 
during and after each shock impulse. Equipment tests of 
this nature are quite expensive and require a great deal of 
preparation. It is desirable for the designer to have some 
idea of what the transient response of the hull and 
equipment will be prior to the actual UNDEX. 

The development of analytical methods and computer codes 
for the analysis of the reaction of submerged structures to 
underwater explosions has enabled the engineer to simulate 
the shock response of submarines and equipment with 
increasing accuracy. Finite-element/finite difference 
methods allow for the analysis of structure and -fluid 
responses, and the work of Geers [Ref. 2] provides a means 
to incorporate fluid-structure interaction effects. Several 
computer codes currently exist which utilize these 
principles to analyze installed equipment response to shock 
waves. One of them is the ELSHOK {^ASTIC SHO CK) code 
developed by Weidlinger Associates, under the sponsorship of 
the Defense Nuclear Agency (DNA) and the Office of Naval 
Research (ONR) . This code was developed to investigate 
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modern submarine underwater explosive shock response in 
conjunction with a testing program using small to large 
scale models and shaped or tapered explosive charges. The 
accuracy of ELSHOK has been validated in several highly 
controlled tests, notably the 1983 low level explosive test 
of an SSN 668 class submarine. The code was made available 
to the Naval Postgraduate School by the DNA with support 
from Weidlinger Associates. It was first used at the school 
by LT Mark Welch, USN , in his study comparing the shock 
response predictions of ELSHOK with those obtained using the 
Dynamic Design Analysis Method (DDAM) [Ref. 3]. 

B. PURPOSE FOR THIS INVESTIGATION 

An explosive detonated underwater exerts a great amount 
of pressure on the surrounding fluid. The water is actually 
compressed by the force of the explosion causing a shock 
wave to form at the point of detonation which propagates in 
a roughly spherical shape at the speed of sound in water. 
When the pressure pulse impinges on a submarine, dynamic 
responses result in the hull from the fluid-structure 
interaction. The submarine experiences translational motion 
away from the point of impact. The velocity of this motion 
is of particular interest to a designer of submarines or 
weapons . 

In an explosion with a side-on aspect with the 
submarine, the shock wave hits the hull in the middle before 
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it hits the ends due to its spherical shape. This causes 
the middle to bow away from the point of impact while the 
ends remain initially fixed. The resulting motion is called 
"whipping" as the hull bends back and forth about the 
middle. With an explosion off the bow or stern, the shock 
wave hits the submarine as a plane wave due to the small 
cross-sectional area at the point of impact. Since the 
force of the shock wave is axisymmetric about the centerline 
of the submarine, no whipping motion is induced in the hull. 
The inertial forces are transmitted down the hull until they 
reach the stern. The stern reacts violently to the combined 
inertial forces resulting from- the explosion as well as the 
reflected pressure pulses from the water displaced by the 
translational motion of the submarine. After the initial 
transient response, the hull settles into an "accordion 
mode" about the middle where the motion of the bow opposes 
that of the stern . 

The purpose of this investigation is to examine the 
shell/substructure interaction when subjected to UNDEX end- 
on loading using the ELSHOK code. The hull used in this 
study is a ring-stiffened cylindrical shell with endplates 
similar in shape to those used as SSTV's. The internal 
equipment or substructure is a diaphragm attached to the 
shell whose thickness is increased in order to study how the 
mass ratio between substructure and shell affects the 
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coupled transient response. The dynamic amplification of 
the response of the substructure at resonance is examined 
for both taper and conventional charges of equivalent 
impulse. 
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II . 



EXPLOSIVE SHOCK RESPONSE USING ELSHOK 



A. GENERAL PRINCIPLES OF OPERATION 

The ELSHOK computer code consists of a family of 
programs developed to calculate the transient response of a 
submerged, ring-stiffened shell of revolution of finite 
length to an underwater shock wave emanating from an 
explosive source located at an arbitrary point away from the 
structure [Ref. 4]. The shell is assumed to be linearly 
elastic, with or without internal substructures, and the 
surrounding fluid is treated as an infinite acoustic medium. 
Component modal analysis is employed in all phases of the 
calculations. The complete structural system is considered 
to consist of the ring-stiffened shell and any attached 
substructures. The- vibration modes of each component are 
calculated separately, and the equations of motion for the 
entire system are obtained by enforcing compatibility of 
deformation at the points of attachment. The free-free 
modes of the empty ring-stiffened shell and the fixed-base 
modes of each individual substructure are coupled through 
the use of dynamic boundary conditions [Refs. 5, 6]. This 
eliminates the need to calculate the combined modes and 
natural frequencies of the entire system as well as the 
requirement for a combined system stiffness matrix. 
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The structure-fluid interaction is approximated in 
ELSHOK using the Doubly Asymptotic Approximation (DAA) 
methods of Geers [Ref. 2], expressed in terms of functions 
which are orthogonal over the wet surface of the submerged 
shell. By matching exact pressure-velocity relations at 
zero and infinite frequencies/ the elements of the matrices 
in the resulting DAA are obtained. In transient problems, 
the DAA yields exact solutions at early and late times 
providing a smooth transition between these two limits. By 
accounting for the effects of the fluid with quantities 
defined solely on the wet surface of the shell, the DAA 
essentially uncouples the fluid field from the structural 
field . 

B. ORGANIZATION AND OPERATION OF ELSHOK 

As Stated previously, the ELSHOK code contains a group 
of computer programs which are utilized in order to obtain 
veloci ty-time histories for various locations in the shell- 
substructure system. The programs, in order of execution, 
are : 

1) B0S0R4 — structural analyzer for shell [Ref. 7] 

2) ACESNID — virtual mass processor 

3) PIFLASH--shell-f luid processor 

4) SAPIV--structural analyzer for substructure [Ref. 8] 

5) PICRUST — substructure processor 

6) USLOB — time integration processor 
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Figure 1. Organization of the ELSHOK Computer Code [Ref. 4] 



7) PUSLOB — plotting processor 
Figure 1 illustrates the general relationship between the 
programs which make up the ELSHOK computer code. As can be 
seen, the operation is divided into four phases. 

1 . Phase I — Shell and Fluid Analysis 

The first step of the ELSHOK calculation is to 
determine the in-vacuo free-free modes and natural frequen- 
cies of the shell. This is accomplished by modeling the 
shell using the B0S0R4 finite difference code. This code 
applies to segmented, ring-stiffened, branched shells of 
revolution having various meridional geometries, wall 
constructions, and ring reinforcements making it well suited 
for the modeling of submarines or SSTV s . For compatibility 
with ELSHOK, a separate B0S0R4 calculation must be performed 

r 

for each circumferential harmonic (N) included in the 
analysis. Examples of different harmonics in the circum- 
ferential distributions include N = 0 (breathing/torsional) 
and N = 1 (pure translation/whipping) as shown in Figure 2. 

The second part of Phase I is the computation of the 
virtual mass array using ACESNID. This provides the 
late-time contribution of the DAA. The virtual mass array 
is determined from the solution, based on simple sources of 
a low-frequency steady-state problem in which normal 
displacements corresponding to surface expansion functions 
are applred to the surface of revolution in the infinite 
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Figure 2. Nodal Patterns for Cylindrical Shells 
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fluid having the same shape and size as the net surface of 
the shell. ACESNID requires only one execution for all 
values of N considered. 

The execution of PIFLASH completes Phase I by 
combining the outputs from B0S0R4 and ACESNID to create a 
"shell-fluid file" which contains all the information 
required to describe the shell during subsequent operations. 
One of the useful features of ELSHOK is the fact that any 
number of substructures can be analyzed for a given shell 
without having to repeat Phase I calculations. 

2 . Phase I I--Substructure Analysis 

The analysis of the substructure is performed using 
the SAPIV finite element code. The SAPIV code contains a 
variety of finite elements (plates, beams, pipes, etc.) 
which can be used in the modeling process. It is a general 
purpose code based upon a lumped mass formulation applying 
to linearly elastic structures. The PICRUST code takes data 
from SAPIV and reorganizes it to facilitate the solution of 
the equations governing the transient response problem under 
study. It is during the execution of PICRUST that the 
connectivity between the substructure and the shell is 
accounted for and influence coefficients for evaluating the 
forces developed at the points of attachment are calculated. 

3 . Phase III — Submerged Shock Response 

In order to prepare a combined shell-substructure 
input for the USLOB time integration processor, the 
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shell-fluid file from Phase I is merged with the 
substructure files from Phase II. USLOB is the portion of 
the ELSHOK code which contains the underwater shock from an 
arbitrary point. It allows pressure-t ime history inputs for 
taper charge modeling or an exponentially decaying pressure 
impulse such as those encountered with conventional charges. 
It is of the form [Ref. 9] : 



K, 



where : 
P(t) 

^1 



P(t) = K^(w^/^/R) ^ exp(-t/eQ) 



incident pressure on the shell .(psi) 

multiplicative constant for incident pressure 

spatial decay constant for incident pressure 

time after arrival of shock wave at point of 
interest (msec) 



( 1 ) 



w 



1/3 1/3 *'’4 

K^w ^ (w ^ /R) = time constant of exponential 



K, 



K, 



weight of spherical charge (lb.) 

K, 

R) 

decay (msec) 

multiplicative constant for time constant 
spatial decay constant for time constant 

R = distance from the explosive to the point of 
interest (ft) 

, K 2 , and are constants which depend on 

explosive type. USLOB employs a modified version of the 
Runge-Kutta integration method to produce veloci ty-time 
histories in tabular form at user specified points on the 
shell and the substructure. 
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4 . Phase IV--Plots of Velocity-Time Histories 
After the velocity-time histories have been 
obtained, PUSLOB is utilized to display the velocity 
responses in a plotted format on a TEKTRONIX graphics 
terminal. These plots are adequate for viewing trends and 
general responses, however for formal display purposes, the 
IBM Versatec plotter yields a superior product. Appendix A 
contains a code which was written to convert the velocity 
punch-card files created by PUSLOB into a data format which 
is acceptable for Versatic plotting using the EASYPLOT 
program. 



24 



III. 



MODEL USED IN THE ANALYSIS 



A. SHELL MODEL 

The shell used in this analysis is a small-scale model 
of a typical SSTV used for conducting underwater shock tests 
on high-impact shipboard machinery. As can be seen in 
Figure 3, it is a high strength steel free-free ring 
stiffened cylinder with aluminum endplates. Appendix B 
includes a sample B0S0R4 input for one of the 
circumferential harmonic distributions. Data can be input 
to B0S0R4 using a formatted input code or by responding 
interactively to a set of prompts. Each segment is modeled 
separately which allows for different materials and physical 
properties to be included in the structure. The small 
stiffeners along the shell are represented by an orthotropic 
approximation which in effect increases the density of the 
shell. The six discrete rings are modeled individually, 
allowing for different size or material. The torsional 
rigidity (GJ) for each ring is found using the following 
relation : 



GJ = EJ/2 (1 + V) 



( 2 ) 



where : 

E = modulus of elasticity (psi) 
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Figure 3. Ring Stiffened Cylindrical Shell Used in Analysis 



4 

J = A /(40 I ) = torsional rigidity constant [Ref. 10] 

V = Poisson' s ratio 

Since the shell is a surface of revolution, the 
properties of each ring cross-sect ion need only be specified 
at one point. The end plates (segments one and three) are 
divided into ten equally spaced nodes, and the cylindrical 
portion of the shell (segment two) is divided into forty- 
five nodes along the longitudinal axis. As with any finite 
difference code, greater accuracy will be obtained in the 
solution with higher numbers of nodes, however, computa- 
tional time will increase. Two additional nodes are 
automatically inserted by BOSOR4 into each segment in order 
to reduce the truncation errors associated with segment 
interfaces and to prevent spurious vibration or buckling 
modes. Through finite difference techniques, BOSOR4 
determines the in-vacuo free-free modes and natural 
frequencies of the shell. Since this study deals strictly 
with end-on loading, the problem is axisymmetric about the 
longitudinal axis. Although only the N = 0 (breathing/ 
torsional) modes are activated, a few N = 1 (translational) 
models must be calculated and retained due to the fact that 
the ELSHOK computer code was developed for side-on loading 
and, therefore, requires a translational input. Since there 
is no translation, the solution is not affected by these 
superfluous modes. BOSOR4 generates a shell file which 
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Figure 4. Typical Shell-Substructure Configuration [Ref. 4] 



gives the mode shapes for each natural frequency. All 
breathing modes have been retained for further study, 
however, the torsional modes have been dropped from the 
analysis since they are not activated by an end-on load. 
Table I contains a listing of modes and natural frequencies 
retained for input into the PIFLASH program. 

Figure 4 gives a representation of a typical shell- 
substructure configuration. The two coordinate systems used 
in an ELSHOK analysis are illustrated. X, Y, and Z refer to 
the global coordinate system of the shell, while x, y, and z 
refer to the local coordinate system of the substructure. 

In both systems, the x-axes run longitudinally down the 
hull. The substructure z-axis is related to the global 
Z-axis by an angle ct^ where a refers to the substructure. 
Points on the wet surface of the shell are located by 
specifying the meridional arc length s, measured along the 
reference surface, and the circumferential coordinate angle 
9. The lower case u, v, and w are used for local shell 
displacements. To eliminate torsional modes from 
consideration, all N = 0 modes in the shell file with 
significant v-d isplacements are discarded. 

Once the shell files have been obtained from BOSOR4, the 
input code for ACESNID can be created. This program 
determines the virtual mass array which produces the late- 
time contribution of the DAA. It considers a cavity in an 



29 



6 

7 

9 

12 

13 

14 

19 

20 

21 

22 

23 

25 

26 

27 

28 

29 

30 

1 

2 



TABLE I 

MODES AND NATURAL FREQUENCIES OF SHELL 

N FREQUENCY (HZ) 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 



1. 259E-03 
2.239E+02 
3.850E+02 
4.761E+02 
5.935E+02 
7.826E+02 
9.369E+02 
1.048E+03 
1. 584E+03 
1.585E+03 
1. 596E+03 
1.598E+03 
1. 656E+03 
1.791E+03 
1. 818E+03 
1.856E+03 
1. 878E+03 
1.890E+03 
1. 898E+03 
1.390E-04 
6. 321E+01 
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infinite acoustic fluid corresponding to the wet surface of 
the shell. A user specified number of surface expansion 
functions having the property of orthogonality over the wet 
surface of the shell are obtained during the calculation of 
the virtual or entrained mass. The cavity is divided into 
bands which describe the behavior of the fluid. A 
sufficient number of bands must be provided to enable the 
fluid response to match the normal motion specified by the 
surface expansion functions. In a sense, they Gan be looked 
upon as nodes in a finite element model. In this analysis, 
the shell is assumed to be totally immersed in salt water 
with a mass density of 9.59684E-05 Ibf-s^/in*^. The 
endplates are divided into fifteen bands each, and the 
cylinder has fifty-one bands. Six surface expansion 
functions have been generated for the endplates with ten 
along the cylinder. 

The shell mode files from B0S0R4 and the virtual mass 
file from ACESNID are combined in the PIFLASH program to 
produce a shell-fluid file for further processing. The 
model of the shell is now complete and can be used for any 
variety of substructures without any recalculation. 

B. SUBSTRUCTURE MODEL 

The idealized internal equipment selected for study in 
this analysis is a high strength steel diaphragm or plate 
located at the discrete ring at frame nine. The SAPIV 
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finite element code is utilized to model the diaphragm. The 
diaphragm has been divided into eight "pie-shaped" wedges 
which are further subdivided into three sections. Twenty- 
four plate elements are used in this model with twenty-five 
nodal points. Since this study is concerned with end-on 
loading, the node at the center of the plate is constrained 
to move only in the x-direction with no rotations. The 
outer nodes (eighteen through twenty-five) are rigidly 
attached to the shell. This is accomplished through a 
modification of the general purpose SAPIV code by the 
developers of ELSHOK. During the parametric study, the 
diaphragm thickness is uniformly varied from .25 inches to 
9.25 inches. Figure 5 is a TEKTRONIX representation of the 
finite element model of the substructure. Appendix B 
includes one of the SAPIV input codes used during the 
parametric study. 
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IV. ANALYSIS 



The purpose of this analysis is to study the shell/ 
substructure interaction after the system is subjected to an 
end-on load while varying the mass and stiffness of the 
internal substructure. Two different types of loading are 
explored: a taper charge and a conventional charge. Taper 

charges are characterized by sustained incident pressure on 
a body over a period of time. They are used in underwater 
shock tests to simulate the type of pressure profile which a 
body would experience from a nuclear detonation. Extensive 
use of scaling is involved in efforts to reproduce the 
physical behavior of an object in response to a nuclear 
charge without generating the extremely high pressures 
associated with actual detonations. Conventional charges 
create a large initial pressure pulse followed by a rapidly 
decaying exponential. ELSHOK uses inputs of charge weight 
and standoff geometry to calculate the transient velocities 
of the shell and substructure in response to these loads. 

A. ANALYSIS PROCEDURE 

Once the shell model has been produced, a series of 
SAPIV runs are conducted for the diaphragms of increasing 
thickness. Each run generates natural frequencies and 
accompanying mode shapes for each model. The PICRUST 
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program reduces the data obtained through SAPIV into a 
format which can be combined with the shell-fluid file for 
processing with the USLOB code. The substructure is 
connected to the shell during this phase of the ELSHOK 
analysis. The diaphragm is attached to shell segment two at 
frame nine (node eleven) . 

The weight, type, and location of the charge is 
specified in the USLOB input code. To set up an end-on load 
referenced to 'the global coordinate system, the charges are 
placed in the negative X-direction at a distance of 70 feet 
(840 inches). This distance is typical of stand-off 
distances used during actual tests. For the taper charge, a 

n Figure 6 is 
rder to make a 
espouses from 
s necessary 
In the field 
to the time- 
e impulse of 
t after its 

(3) 

Iculating the 
6. As was 



typical incident pressure- time history shown i 
used to simulate the shock wave loading. In o 
meaningful comparison between the structural r 
a taper charge and a conventional charge, it i 
that they have equivalent shock wave impulses, 
of underwater shock analysis, "impulse" refers 
integral of the pressure profile [Ref. 9]. Th 
unit area of the shock wave front up to a time 
arrival is given by: 

I(t) = / P(t)dt 

The impulse of the taper charge is found by ca 
area under the pressure-time history in Figure 
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PRESSURE-TIME HISTORY 



TAPER CHARGE 




Figure 6« Incident Pressure-Time History — Taper Charge 

PRESSURE-TIME HISTORY 



CONVENTIONAL CHARGE 




Figure 7. Incident Pressure-Time History — Conventional 
Charge 
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stated earlier, a conventional charge follows a pressure- 
time history given by equation (1) : 



P(t) = ^ exp(-t/e^) 



( 1 ) 



where : 




time constant of exponential 
decay (msec) 



HBX-1 is a popular explosive used in underwater shock tests, 
so it was decided to use it in this study for the 
conventional charge. K 2 / , and are constants which 

depend on the type of explosive used. The values used in 
this study are [Ref. 11]: 

= 3.8354314E+05 

K 2 = 1.144 

= 3.03131E-05 

= -.247 

By matching the impulses between the taper charge and the 
conventional charge, a charge weight of 352 lbs of HBX-1 is 
found to yield an equivalent impulse. The peak pressure 
generated by the charge is 1620 psi which creates the 
incident pressure-time history shown in Figure 7. The 
conventional charge expends most of its impulse in the first 
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millisecond, while the taper charge takes over four milli- 
seconds to expend an equivalent amount of impulse. Because 
of the sharp rise in pressure experienced from a conven- 
tional charge, a higher transient response is developed in 
the shell than with a taper charge of equivalent impulse. 

After selecting suitable charges, the time step incre- 
ment and integration limits are specified in the USLOB input 
code. Through a trial and error process, enough time steps 
are chosen to identify the significant interactions between 
the shell and substructure. It was found that 1600 time 
steps for an 80 msec time period gave an adequate illustra- 
tion of the system response. The USLOB program enables the 
user. to examine the velocity-time history response of 
various nodes in tabular form. This gives some indication 
of whether the run was successful, however, the plots from 
PUSLOB are required before response characteristics can be 
ident i f ied . 

To obtain a baseline velocity-time history, a complete 
analysis is performed on the submerged empty stiffened 
shell. From the velocity-time history plots, the dominant 
frequency from the computed shock response is observed. 
Provided the mass of the diaphragm remained small with 
respect to the shell, it was felt that this provided a good 
estimation of the excitation frequency for the diaphragm. 

The dominant frequency is 187.80 Hz for the taper charge and 
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195.44 Hz for the conventional charge. A series of SAPIV 
calculations are conducted on diaphragms of varying 
thickness in order to identify the fundamental frequency of 
each substructure. These values are given in Table II. An 
analysis is performed on the shell/substructure system for 
each diaphragm thickness to observe the coupled velocity- 
time history response as the substructure mass grows in 
relation to that of the shell. 
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TABLE II 



FUNDAMENTAL FREQUENCIES 

Thickness (In) 

.25 

.50 

.75 

1.00 

1.05 

1.06 
1.07 
1.10 
1.15 

1.25 

2.25 

3.25 

4.25 

5.25 

6.25 

7.25 

8.25 

9.25 



OF DIAPHRAGM SUBSTRUCTURE 
Frequency (Hz) 



43 , 


. 66 


87 , 


.31 


130 , 


.97 


174 , 


.59 


183 , 


. 35 


185 , 


.10 


186 . 


. 85 


192 . 


.10 


200 . 


, 85 


218 . 


.20 


392 . 


, 95 


567 . 


,55 


742 . 


, 14 


916 . 


,73 


1091 . 


33 


1265 . 


,92 


1440 . 


51 


1615 . 


,42 
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V 



RESULTS 



The velocity time-history response begins at time zero 
when the incident pressure wave from the underwater 
explosion strikes the shell. With end-on loading, the 
forward endplate receives the initial contact with velocity 
indicated along the longitudinal axis. The empty shell 
response is analyzed first, followed by a discussion of the 
coupled effects when the diaphragm achieves resonance. 

A. VELOCITY-TIME HISTORY RESPONSE — EMPTY SHELL 

The empty shell velocity-time history responses for the 
conventional and taper charges are contained in Figures 8 
through 20. It can be seen that the conventional charge 
causes a much greater transient response in the shell. In 
both cases, the velocity is rapidly damped by interaction 
between the shell and the infinite fluid medium. It is 
interesting to observe the differences in response between 
the various locations on the shell. The forward endplate is 
perpendicular to the shock wave, receiving the full force 
from the incident pressure pulse. Inertial forces are 
generated by the endplate causing severe local deformation 
in the nearest nodes on the cylinder due to the rigid 
coupling between segments. The motion propagates down the 
shell along the longitudinal axis until it reaches the aft 
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VELOCITY RESPONSE - EMPTY SHELL 

SEGMENT ONE - NODE ONE 
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Figure 8. Velocity-Time History Response (Empty Shell) — FWD Endplate Center 



VELOCITY RESPONSE - EMPTY SHEU. 

SEGMENT TKO - NODE TWO 
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Figure 9. Velocity-Time History Response (Empty Shell) — Frame Zero 



VELOCITY RESPONSE - EMPTY SHELL 

SEGMENT TWO - NODE SEVEN 
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Figure 10. Velocity-Time History Response (Empty Shell) — Frame Five 



VELOCITY RESPONSE - EMPTY SHELL 

SEGMENT TWO - NODE ELEVEN 
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Figure 11. Velocity-Time History Response (Empty Shell) — Frame Nine 



VELOCITY RESPONSE - EMPTY SHELL 

SEGMENT TWO - NODE FIFTEEN 
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Figure 12. Velocity-Time History Response (Empty Shell) — Frame Thirteen 



VELOCITY RESPONSE - EMPTY SHELL 

SEGMENT TWO - NODE TWENTY 




(D3S/Xi) A1I30T3A 
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Figure 13. Velocity-Time History Response (Empty She 11 ) --Frame Eighteen 



VELOCITY RESPONSE - EMPTY SHELL 

SEGMENT TWO - NODE TWENTY FOUR 




(D2S/1J) AXD013A 
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Figure 14. Velocity-Time History Response (Empty Shell) — Frame Twenty-Three 



VELOCITY RESPONSE - EMPTY SHELL 

SEGMENT TWO - TWENTY EIGHT 
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Figure 15. Velocity-Time History Response (Empty Shell) — Frame Tv/enty-Seven 



VELOCITY RESPONSE - EMPTY SHELL 

SEGMENT TWO - NODE THIRTY THREE 




50 



Figure 16. Velocity-Time History Response (Empty Shell) — Frame Thirty-Two 



VELOCITY RESPONSE - EMPTY SHELL 

SEGMENT TWO - NODE THIRTY SEVEN 
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Figure 17. Velocity-Time History Response (Empty Shell) --Frame Thirty-Six 



VELOCITY RESPONSE - EMPTY SHELL 

SEGMENT TWO - NODE FORTY ONE 
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Figure 18. Velocity-Time History Response (Empty Shell) — Frame Forty 



VELOCITY RESPONSE - EMPTY SHELL 

SEGMENT TWO - NODE FORTY SIX 




(D3S/ii) AlIOOiaA 
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Figure 19. Velocity-Time History Response (Empty Shell) — Frame Forty-Five 



VELOCITY RESPONSE - EMPTY SHELL 

SEGMENT THREE - NODE TWELVE 




(oas/xj) AXDoiaA 
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Figure 20. Velocity-Time History Response (Empty Shell) — AFT Endplate Center 



endplate over a millisecond after initial shock wave 
contact. At first the aft endplate moves in the negative 
X-direction as the "ripple" proceeds down the structure. 

Then it experiences a violent transient response as it is 
affected not only by the inertial forces from the shell, but 
by the reflected pressure pulses from the water it 
displaces. After the initial shock, the shell settles into 
an "accordion" effect, expanding and contracting as the 
forward and aft endplates move in opposite directions. 

B. VELOCITY-TIME HISTORY RESPONSE — COUPLED SHELL/ 

SUBSTRUCTURE SYSTEM 

From the results of the analyses performed on the coupled 
shell/substructure system. Tables III and IV were created in 
order to identify when the diaphragm achieved resonance. 
is the frequency ratio between the fundamental frequency of 
the simply-supported diaphragm (Q) and the dominant frequency 
observed in the veloci ty-t ime history response of the sub- 
merged empty stiffened shell (ii) . M1/M2 is the mass ratio 

between the substructure (M^^) and the shell with endplates 
(M2). Figure 21 is a plot of the velocity of the diaphragm 
center versus the frequency ratio Si/S. it reveals that in 
both the taper and conventional charge responses, the 
diaphragm achieves its maximum velocity or resonance when the 
fundamental frequency of the substructure is nearly equal to 
the excitation frequency from the shell. Resonance occurs 
when the mass ratio is approximately .08. Since 
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TABLE III 



MAXIMUM VELOCITY RESPONSE OF DIAPHRAGM SUBSTRUCTURE 
WITH VARYING THICKNESS (TAPER CHARGE) 



Thickness (In) 




M1/M2 


Velocity ( 


.25 


.232 


.018 


22.19 


.50 


.465 


.036 


16.64 


.75 


.762 


.055 


20.79 


1.00 


.930 


.073 


34.81 


1.05 


.976 


.076 


35.65 


1.06 


.986 


.077 


35.78 


1.07 


.995 


.078 


35.65 


o 
1 — 1 

r — 1 


1.023 


.080 


33.58 


1.15' 


1.069 


. 084 


30.36 


1.25 


1.162 


.091 


26.24 


2.25 


2.092 


. 164 


11.61 


3.25 


3.022 


.237 


10.60 


4.25 


3.952 


.310 


9.04 


5.25 


4.881 


.383 


8.70 


6.25 


5.811 


.456 


8.44 


7.25 


6.741 


.528 


8.15 


8.25 


7.670 


.601 


7.87 


9.25 


8.602 


.674 


7.62 
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TABLE IV 



MAXIMUM VELOCITY RESPONSE OF DIAPHRAGM SUBSTRUCTURE 
WITH VARYING THICKNESS (CONVENTIONAL CHARGE) 



Thickness (In) 




M1/M2 


Velocity ( 


.25 


.223 


. .018 


40.42 


.50 


.447 


.036 


28.26 


.75 


.670 


.055 


45.29 


1.00 


.893 


.073 


68.73 


1.05 


.938 


.076 


73.84 


1.06 


.947 


.077 


' 74.71 


1.07 


.956 


.078 


74.22 


1.10 


.983 


.080 


69.69 


1.15 


1.028 


.084 


62.90 


1.25 . 


1.116 


.091 


62.75 


2.25 


2.011 


.164 


21.77 


3.25 


2.904 


.237 


17.67 


4.25 


3.797 


.310 


15.51 


5.25 


4.691 


• .383 


13.69 


6.25 


5.584 


.456 


12.01 


7.25 


6.477 


.528 


11.37 


8.25 


7.371 


.601 


11.00 


9.25 


8.265 


.674 


10.69 
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DIAPH MAX VEL VS NORMALIZED FREQ 




(oas/u) aiND RovaHdvia aaA xvn 
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Figure 21. Diaphragm Maximum Velocity vs. Normalized Frequency 



SHELL RESPONSE 
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Figure 22. Two Degrees of Freedom System 
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(4) 



SJ = JJ = V Kl/Ml = ^K2/M2' 

where: K1 = overall stiffness of the substructure 
K2 = overall stiffness of the shell 

it can be said that: 

Kl/Ml = K2/M2 = K1/.08 M2 

or 

K1/K2 = .08 

The coupled shell and substructure can then be modeled as a 
two degrees of freedom system as shown in Figure 22 where Ml 
and Kl represent the substructure and M2 and K2 represent 
the shell. The large mass (M2) has an uncoupled natural 
frequency (£ 2 ) which is equal to aJ K2/M2'. This is 
equivalent to the dominant frequency (j^) of the empty shell. 
Since this is a two degrees of freedom system, it will 
possess two natural frequencies for the coupled system. It 
has been found in coupled systems that for decreasing mass 
ratio (M1/M2) and for increasing stiffness ratio (K1/K2) , 
the first and second modal frequencies of a two degrees of 
freedom system approach the decoupled frequencies f^ and 
of two single degrees of freedom systems, respectively [Ref. 
12]. The first natural frequency is bounded by 
is bounded by the second natural frequency. From Figure 23, 
it is seen that with M1/M2 = K1/K2 = .08, the dominant 
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frequency of the shell is approximately equal to the first 
modal frequency of the system. The result is that the 

shell/substructure system achieves resonance when excited at 
the dominant frequency of the empty shell. 

With the diaphragm at resonance, Figures 24 and 25 
illustrate that the transient response of the diaphragm is 
much more pronounced than that of the shell due to the large 
differences between the two masses as well as the soft 
spring between them. Initially the diaphragm is drawn in 
the negative X-direction by inertial forces as the shell is 
displaced by the shock wave. With the increased motion of 
the shell at frame nine, the diaphragm receives more and 
more energy. At approximately t = 20 msec, the diaphragm is 
at its maximum velocity while the shell is relatively calm. 
At this point, the diaphragm begins to transfer energy back 

p 

to the shell. Since kinetic energy is equal to 1/2 MV”, the 
diaphragm is a small mass with high velocity while the much 
larger mass of the shell moves at a lower velocity. The 
energy of the system is transferred between the shell and 
substructure as the overall motion is damped by interaction 
between the shell and the fluid medium. The transfer of 
energy between the shell and the diaphragm creates an effect 
known as the "beating" phenomenon. This effect is seen in 
the velocity time history responses for the case of 
diaphragm resonance in Figures 26 through 39. The effect is 
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VELOCITY RESPONSE - T=1.06 IN. 

TAPER CHARGE 




(Das/ii) A1ID013A 
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Figure 24. Energy Exchange Between Shell and Diaphragm (Taper Charge) 



VELOCITY RESPONSE - T=1.06 IN. 

CONVENTIONAL CHARGE 




(oas/li) AlIDOiaA 
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Figure 25. Energy ' Exchange Between Shell and Diaphragm (Conventional Charge) 



VELOCITY RESPONSE - T=1.06 IN. 

SEGMENT ONE - NODE ONE 
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Figure 26. Velocity-Time History Response (Resonant Case) --FWD Endplate Center 



VELOCITY RESPONSE - 1^1.06 IN. 

SEGMENT TWO - NODE TWO 
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Figure 27. Velocity-Time History Response (Resonant Case) --Frame Zero 



VELOCITY RESPONSE - 'Ii=1.06 IN. 

SEGMENT TWO - NODE SEVEN 
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Figure 28. Velocity-Time History Response (Resonant Case) --Frame Five 



VELOCITY RESPONSE - T=1.06 IN. 

SKGMENT TiyO - NODE ELEVEN 
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Figure 29. Velocity-Time History Response (Resonant Case) — Frame Nine 



VELOCITY RESPONSE - T=1.06 IN. 

DIAPHRACM CENTER 
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Figure 30. Velocity-Time History Response (Resonant Case) — Diaphragm Center 



VELOCITY RESPONSE - 'P=1.06 IN. 

SEGMENT TWO - NODE FIFTEEN 




70 



Figure 31. Velocity-Time History Response (Resonant Case) — Frame Thirteen 



VELOCITY RESI’ONSE - '1^1.06 IN. 

SEGMENT TWO - NODE TWENTY 
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Figure 32. Velocity-Time History Response (Resonant Case) — Frame Eighteen 



VELOCITY RESPONSE - T=1.06 IN. 

SEGMENT TWO - NODE TWENTY FOUR 
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Figure 33. Velocity-Time History Response (Resonant Case) — Frame Twenty-Three 



VELOCITY RESPONSE - T=1.06 IN. 

SEGMENT TWO - NODE TWENTY EIGHT 
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Figure 34 . Velocity-Time History Response (Resonant Case) Frame Twenty Seven 



VELOCITY RESPONSE - Ti=1.06 IN. 

SEGMENT TWO - NODE THIRTY THREE 




74 



Figure 35. Velocity-Time History Response (Resonant Case) — Frame Thirty-Two 



VELOCITY RESPONSE - 'V=L06 IN. 

SEGMENT TWO - NODE THIRTY SEVEN 
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Figure 36. Velocity-Time History Response (Resonant Case) — Frame Thirty-Six 



VELOCITY RESPONSE - T=1.06 IN. 

SEGMENT TWO - NODE FORTY ONE 
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Figure 37. Velocity-Time History Response (Resonant Case) — Frame Forty 



VELOCITY RESPONSE - 'Ii=1.06 IN. 

SEGMENT TWO - NODE FORTY SIX 




(33S/li) AXI3013A 
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Figure 38. Velocity -Time History Response (Resonant Case) --Frame Forty-Five 



VELOCITY RESPONSE - 'r=1.06 IN. 

SEGMENT THREE - NODE TWELVE 
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Figure 39. Velocity-Time History Response (Resonant Case) — AFT Endplate Center 



quite pronounced at the endplates, because they act as "hard 
spots" receiving the energy transmitted through the shell. 
The cylindrical portion of the shell acts as a conduit for 
the energy transfer. Again, the nodes on the shell nearest 
the endplate experience some local deformation due to the 
rigid coupling between segments. 

The empty shell velocity-time history response is 
compared with the resonant diaphragm case in Figures 40 
through 45. Of particular interest is the fact that the 
velocity profiles are virtually identical for the first 
fifteen milliseconds. During this period of time, the shell 
is receiving energy from the shock wave and transferring a 
portion of it to the diaphragm. Because of the differences 
in mass between the shell and substructure, not much energy 
is lost from the shell. The peak velocities do not vary by 
more than five feet per second during this time period. The 
difference arises at t = 35 msec when the shell receives 
energy from the resonating plate. The decaying velocity 
response is disturbed during the exchange of energy which- is 
observed in the beating effect. 

As the thickness of the diaphragm is increased beyond 
the resonant condition, the interactions between the shell 
and substructure are reduced dramatically. The shell/ 
substructure velocity-time history remains within 10% of the 
empty shell response until the diaphragm thickness is 2.25 
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EMPTY SHELL VS. T^LOe IN. 

SEGMENT ONE - NODE ONE 
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Figure 40. Comparison of Empty vs. Resonant Shell Response (Taper Charge) --Forv/ard 
Endplate Center 



EMPTY SHELL VS. T=1.06 IN. 

SEGMENT TWO - NODE ELEVEN 
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Figure 41. Comparison of Empty vs. Resonant Shell Response (Taper Charge) --Frame 



EMPTY SHELL VS. Ti^LOe IN. 

SEGMENT THREE - NODE TWELVE 




(oas/ii) AUDoaaA 
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Figure 42. Comparison of Empty vs. Resonant Shell Response (Taper Charge) — AFT 
Endplate Center 



EMPTY SHELL VS. T=1.06 IN 

SEGMENT ONE - NODE ONE 



I 




(Das/Xi) AlIDOiaA 
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Figure 43. Comparison of Empty vs. Resonant Shell Response (Conventional Charge) - 
Forward Endplate Center 



EMPTY SHELL VS. T=1.06 IN. 

SEGMENT TWO - NODE ELEVEN 
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Figure 44. Comparison of Empty vs. Resonant Shell Response (Conventional Charge)- 



EMPTY SHELL VS. T=1.06 IN. 

SEGMENT THREE - NODE TWELVE 
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Figure 45. Comparison of Empty vs. Resonant Shell Response (Conventional Charge) - 
AFT Endplate Center 



inches thick or M1/M2 = .164. By the time it is 9.25 inches 
thick, it represents over half the mass of the shell, and it 
acts as if it was an integral part of the shell slaved to 
the motions of the forward endplate. 
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VI . 



CONCLUSIONS 



Upon inspection of the velocity-time history responses 
at various locations along the shell, one can see why the 
potential for damage is so great from an underwater 
explosion off the bow or stern of a submarine. While the 
cylindrical portion for the most part serves to transmit 
inertial forces with little velocity response, the 
extremities of the hull experience violent transient 
responses which may be passed to attached internal equipment 
located in those areas. The inertial influence of the mass 
of the endplates on neighboring nodes of the shell is 
apparent in both the empty shell and resonant cases. The 
response from a conventional charge is seen to be much more 
pronounced than that of a taper charge of equivalent impulse 
due to the rapid expenditure of energy against the forward 
endplate . 

The results obtained from the resonant diaphragm case 
support the theory that with decreasing mass ratio between 
substructure and shell, the first and second modal 
frequencies of the two degrees of freedom system approach 
the decoupled frequencies of the substructure and shell in 
two single degree of freedom systems. This allows the 
motion of hull to drive the diaphragm into resonance when 
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the dominant frequency of the shell motion is roughly equal 
to the fundamental frequency of the substructure. It is 
possible that this could occur during an underwater shock 
test of a small piece of equipment in a large SSTV. This 
might result in the failure of the equipment in a scenario 
which would not occur in a real submarine with its myriad of 
internal masses and substructures. 

The large difference in mass between the substructure 
and the shell results in a tremendous dynamic amplification 
in the response of the diaphragm at resonance. As energy is 
transferred between the two bodies, the large mass of the 
shell causes its velocity to be quite small compared to that 
of the diaphragm. As is seen in Figures 40 through 45, the 
initial resonant case velocity response is virtually 
unchanged from the empty shell response. If the mass of the 
internal equipment was sufficiently small in comparison to 
an SSTV used in a series of shock tests, one ELSHOK calcu- 
lation could be performed on the empty shell to generate a 
set of velocity response curves. The veloci ty-time history 
response at the attachment point could then be input into 
the SAPIV code for the finite element model of any number of 
small pieces of equipment, and the substructure transient 
response could be obtained without having to run a separate 
ELSHOK calculation for each shell/substructure system. This 
would result in considerable savings in computer time. 
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Although ELSHOK is a powerful tool which can be used to 
accurately predict the results of underwater explosion tests 
on equipment, some shortfalls exist in the program which 
limit its use in design applications. At the present time, 
there is no direct method to take the effects of the surface 
cut-off wave or bottom reflection into account. The surface 
cut-off wave forms in shallow underwater explosions when the 
compressive shock wave reflects off of the surface of the 
water causing a tensile wave to propagate. Bottom 
reflections form in deep explosions as the shock wave is 
reflected off of the floor of the ocean. ELSHOK assumes 
that the explosion is occurring in an infinite acoustic 
medium with a single shock wave. Another weakness in the 
program is the fact that it cannot process oblique 
explosions; charges must be placed in a pure end-on or 
side-on configuration. Shock tests conducted in accordance 
with MIL-S-901D are either end-on or side-on loa<Jing, 
however, so this is not a serious drawback in predicting the 
transient motion of equipment in SSTV's. 

It is recommended that further study be conducted in 
this area using more refined shell models resembling full- 
scale submarines. Once a shell model is produced for a 
given submarine class, researchers can create different 
finite element models of actual equipment to observe the 
interactions between the submarine hull and equipment of 
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varying mass and shape when subjected to end-on or side-on 
loads. It is felt that ELSHOK is an effective method to 
give an indication of the response of equipment in 
underwater explosion tests without needless destruction from 
unexpected failures. Its continued use prior to UNDEX tests 
will lead to savings in both time and effort. 
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APPENDIX A 



PROGRAM USED TO CONVERT PUSLOB PUNCH CARD FILES 
TO FORMAT ACCEPTANCE BY EASYPLOT 

Along with the TEKTRONIX plots of the veloc i ty-t ime 
histories, PUSLOB generates velocity punch card files for 
each requested node. The punch card file consists of a 
block of velocities in exponential format the size of which 
depends upon the number of time points requested . The 
EASYPLOT program developed for the IBM 3033 by Mr. John 
Mainwaring is a quick interactive method to produce graphs 
using DISSPLA without having to create a plotting program. 

It will accept up to 100 data points in tabular form 
provided it is in a column format for X-Y input. The 
PLOTCONVl code was created to accept data from tv/o separate 
punch card files for taper and conventional charge 
comparisons and combine them into one data file for each 
node for input into EASYPLOT. 

In order to prevent data overflow, the time step must be 
sized in USLOB and PUSLOB to generate no more than 100 
velocity data points. Both the PLOTCONVl VAX 11/780 VMS 
computer command file and the FORTRAN code are included in 
this appendix. The command file must be modified prior to 
each run to specify the punch card files being read and the 
titles of the new data files. 
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$! PLOTCONVl.COM - - COMMAND FILE FOR EXECUTION OF PLOTCONVl 
$! 

$ ASSIGN PUSL48.PUN FOROll 
$ ASSIGN PUSL49.PUN FOR012 
$ ASSIGN MAXSHL1.DAT FOR021 
$ ASSIGN MAXSHL2.DAT FOR022 
$ ASSIGN MAXSHL3.DAT FOR023 
$ ASSIGN MAXSHL4.DAT FOR024 
$ ASSIGN MAXSHL5.DAT FOR025 
$ ASSIGN MAXSHL6.DAT FOR026 
$ ASSIGN MAXSHL7.DAT FOR027 
$ ASSIGN MAXSHL8.DAT FOR028 
$ ASSIGN MAXSHL9.DAT FOR029 
$ ASSIGN MAXSHL10.DAT FOR030 • 

$ ASSIGN MAXSHL11.DAT FOR031 
$ ASSIGN MAXSHL12.DAT FOR032 
$ ASSIGN MAXSHL13.DAT FOR033 
$ ASSIGN MAXSHL14.DAT FOR034 

$ RUN PLOTCONVl 



PROGRAM PLOTCONVl 

THE PURPOSE OF THIS PROGRAM IS TO READ TWO PUNCH-CARD 
FILES FROM PUSLOB AND CONVERT THEM TO A FORMAT WHICH CAN 
BE ACCEPTED BY EASYPLOT. 

DIMENSION VEL(100,14) ,VEL0(1 00,3), VELK 100, 14) , LABEL (14), TAG(14), 
%LTITLE(14) 

CHARACTER3€10 LABEL 

CHARACTER5C40 LTITLE 

COMMON/CSOLVE/DTRECS 

NCURV=14 

DO 10 K=1,NCURV 

READ(11,999) LTITLEC K) , NRECS , DTRECS 
READ(11,998) LABEL(K), TAG(K) 

READ(11,997) (VEL(J,K), J=l, NRECS) 

READ(12,999) LTITLE(K) , NRECS, DTRECS 
READ(12,998) LABEL (K) , TAGCK) 

READ(12,997) (VEL1(J,K), J=l, NRECS) 

10 CONTINUE 

999 FORMAT(A40,I5,1PE11.4) 

998 FORMAT(A10,F10.5) 

997 F0RMAT(1P7E11.4) 

DO 20 L=1,NCURV 

T = 0.0 

LL=L+20 

DO 30 M=l, NRECS 
VELO(M,l)=T 
DTRECS1=1000 . OXDTRECS 
T=T+DTRECS1 

VEL0CM,2)=VEL(M,L)/12. 0 
VELO(M,3)=VEL1(M,L)/12.0 
30 CONTINUE 

DO 40 KK=1, NRECS 

WRITE(LL,996) ( VELOC KK, J J ) , J J=1 , 3) 

40 CONTINUE 

WRITE(LL,995) LTITLEC L) , NRECS, DTRECS 
WRITE(LL,994) LABEL ( L ), TAG( L) 

20 CONTINUE 
996 F0RMAT(3F15.5) 

995 F0RMAT(A40,I5,1PE11 .4) 

994 FORMAT(A10,F10 .5) 

STOP 
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APPENDIX B 



TYPICAL ELSHOK INPUT CODES — SHELL AND SUBSTRUCTURE 

The input codes which follow are for the resonant 
diaphragm case for both taper and conventional charges. The 
only differences between the two is in the USLOB code where 
the type of charge is specified. All the input codes are 
explained in detail in the ELSHOK users manual [Ref. 4]. 

1. BOSOR4 INPUT DATA 

The B0S0R4 input code is created by working inter- 
actively with the computer. The prompt following the $ 
symbol appears, and the user supplies the required 
information. The three segments (two endplates and 
cylinder) are modeled separately, and then they are 
"connected" in the global data section. A separate input 
code is required for each circumferential harmonic included 
in the calculation. B0S0R4 output provides the user with 
the in-vacuo free-free modes and natural frequencies of the 
shell . 
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FILE: BOSSO DATA A1 



0S0R4 ATTEMPT ONE N=0 



10 

3 

1 

0.0 

0.0 

16.81250 

0.0 

0 

3 

2.500000 
N 

0 

O.OOOOOOOE+00 

0 

0 

2 

0.1080000E+08 
0 .3200000 
0.2535000E-03 
O.OOOOOQOE+00 
0 
1 
N 
N 
H 
H 

45 

3 

1 

16.81250 

O.OOOOOOQE+00 

16.81250 

258.1250 

0 

3 

0.1875000 

N 

6 

2 

3.500000 
53.12500 
103.7500 
154.3750 
205.0000 
254.6250 

1 

1 

1 

1 

1 

1 

0.3000000E+08 



INDIC = ANALYSIS TYPE INDICATOR 

NPRT = OUTPUT OPTIONS (1=MINIMUM, 2=MEDIUM, 3=MAXIMUM) 

ISTRES= OUTPUT CONTROL ( 0=RESULTANTS » 1=SIGMA. 2=EPSIL0N) 

NSEG = NUMBER OF SHELL SEGMENTS (LESS THAN 25) 

SEGMENT NUMBER 11111111 
NMESH = NUMBER OF NODE POINTS (5 = MIN.; 98 = MAX . ) ( 1) 
NTYPEH= CONTROL INTEGER (1 OR 2 OR 3) FOR NODAL POINT SPACING 
NSHAPE= INDICATOR (1,2 OR 4) FOR GEOMETRY OF MERIDIAN 
R1 = RADIUS AT BEGINNING OF SEGMENT (SEE P. 66) 

21 = AXIAL COORDINATE AT BEGINNING OF SEGMENT 

R2 = RADIUS AT END OF SEGMENT 

22 = AXIAL COORDINATE AT END OF SEGMENT 

IMP = INDICATOR FOR IMPERFECTION (Q=NONE, l=SOME) 

NTYPEZ= CONTROL (1 OR 3) FOR REFERENCE SURFACE LOCATION 
2VAL = DISTANCE FROM LEFTMOST SURF. TO REFERENCE SURF. 

DO YOU WANT TO PRINT OUT R(S), R»(S), ETC. FOR THIS SEGMENT? 
NRINGS- NUMBER (MAX=20) OF DISCRETE RINGS IN THIS SEGMENT 
K=ELASTIC FOUNDATION MODULUS (E.G. LB/IN*x3)IN THIS SEG. 
LINTYP= INDICATOR (0, 1, 2 OR 3) FOR TYPE OF LINE LOADS 
NLTYPE=CONTROL (0,1, 2, 3) FOR TYPE OF SURFACE LOADING 
NWALL=INDEX (1, 2, 4, 5, 6, 7, 3) FOR WALL CONSTRUCTION 
E = YOUNG'S MODULUS FOR SKIN 

U = POISSON'S RATIO FOR SKIN 

SM =MASS DENSITY OF SKIN (E.G. ALUM. = .00025 LB-SECX*2/IN3<5<4 ) 
ALPHA = COEFFICIENT OF THERMAL EXPANSION 

ANRS = CONTROL (0 OR 1) FOR ADDITION OF SMEARED STIFFENERS 
SUR = CONTROL FOR THICKNESS INPUT (0 OR 1 OR -1) 

DO YOU WANT TO PRINT OUT THE C(I,J) AT MERIDIONAL STATIONS? 

DO YOU WANT TO PRINT OUT DISTRIBUTED LOADS ALONG MERIDIAN? 



SEGMENT 
NMESH = 
NTYPEH= 
NSHAPE= 
R1 
21 
R2 
22 
IMP 
NTYPEZ= 
2VAL = 
DO YOU 
NRINGS= 
NTYPE = 
2(1) = 
2 ( 1 ) = 
2 ( 1 ) = 
2 ( 1 ) = 
2 ( 1 ) = 
2 ( 1 ) = 
NTYPER= 
NTYPER= 
NTYPER= 
NTYPER= 
NTYPER= 
NTYPER= 
E 



NUMBER 22222222 
NUMBER OF NODE POINTS (5 = MIN.; 98 = MAX.)( 2) 

CONTROL INTEGER (1 OR 2 OR 3) FOR NODAL POINT SPACING 
INDICATOR (1,2 OR 4) FOR GEOMETRY OF MERIDIAN 
= RADIUS AT BEGINNING OF SEGMENT (SEE P. 66) 

= AXIAL COORDINATE AT BEGINNING OF SEGMENT 
= RADIUS AT END OF SEGMENT 
= AXIAL COORDINATE AT END OF SEGMENT 
INDICATOR FOR IMPERFECTION (0=NONE, l=SOME) 

CONTROL (1 OR 3) FOR REFERENCE SURFACE LOCATION 
DISTANCE FROM LEFTMOST SURF. TO REFERENCE SURF. 

WANT TO PRINT OUT R(S), R'(S), ETC. FOR THIS SEGMENT? 
NUMBER (MAX=20) OF DISCRETE RINGS IN THIS SEGMENT 
CONTROL FOR IDENTIFICATION OF RING LOCATION (2=2, 3=R) 
AXIAL COORDINATE OF ITH RING, 2( 1) 

AXIAL COORDINATE OF ITH RING, 2( 2) 

AXIAL COORDINATE OF ITH RING, 2( 3) 

AXIAL COORDINATE OF ITH RING, 2( 4) 

AXIAL COORDINATE OF ITH RING, 2( 5) 

AXIAL COORDINATE OF ITH RING, 2( 6) 



TYPE (0 OR 
TYPE (0 OR 
TYPE (0 OR 



OR 

OR 

OR 

OR 

OR 

OR 



OR 5) OF DISCRETE RING NO . ( 1) 

OR 5) OF DISCRETE RING NO . ( 2) 

OR 5) OF DISCRETE RING NO . ( 

OR 5) OF DISCRETE RING NO . ( 

OR 5) OF DISCRETE RING NO . ( 5) 

OR 5) OF DISCRETE RING NO . ( 6) 



3) 

4) 



7.125000 


$ 


A 


7.535500 


$ 


lY 


2.375000 


$ 


IX 


0. OOOOOOOE+00 


$ 


IXY 


-1.781250 


$ 


El 


0 . OOOOOOOE+00 


$ 


E2 


0.7501E+08 


$ 


GJ 


0.7330000E-03 


$ 


RM 


0.3000000E+08 


$ 


E 


3.250000 


$ 


A 


11.69800 


$ 


lY 


2.750000 


$ 


IX 


O.OOOOOOOE+00 


$ 


IXY 



TYPE (0 OR 1 OR 
TYPE (0 OR 1 OR 
TYPE (0 OR 1 OR 
YOUNG'S MODULUS OF RING( 1) 

CROSS SECTION AREA OF RING( 1) 

MOMENT OF INERTIA ABOUT Y-AXIS (SEE FIG. ON P.70)( 1) 
MOMENT OF INERTIA ABOUT X-AXIS( 1) 

PRODUCT OF INERTIA( 1) 

RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 70)( 1) 
AXIAL COMPONENT OF RING ECCENTRICITY( 1) 

TORSIONAL RIGIDITY( 1) 

RING MATERIAL DENSITY (E.G. ALUMINUM^ . 00 02535 ) ( 1) 
YOUNG'S MODULUS OF RINGC 2) 

CROSS SECTION AREA OF RING( 2) 

MOMENT OF INERTIA ABOUT Y-AXIS (SEE FIG. ON P.70)( 2) 
MOMENT OF INERTIA ABOUT X-AXIS( 2) 

PRODUCT OF INERTIA( 2) 
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FILE: BOSSO DATA A1 



-2.062500 
O.OOOOOOOE+00 
0 .92A9000E+08 
0.7330000E-03 
0.30Q0000E+08 
' 8.250000 
11.69800 
2.750000 
O.OOOOOOOE+00 
-2.062500 
O.OOOOOOOE+00 
0.92A9000E+08 
0.7330000E-03 
0 .3000000E+08 

8.250000 
11.69800 
2.750000 

O.OOOOOOOE+00 

-2.062500 

O.OOOOOOOE+00 

0.92A9000E+08 

0.7330000E-03 

0.3000000E+08 

3.250000 
11.69800 
2.750000 

O.OOOOOOOE+00 
-2.062500 
O.OOOOOOOE+00 
0.92A9000E+08 
0.7330000E-03 
0.3000000E+08 
7.125000 
7.535500 
2.375000 
0. OOOOOOOE+00 
-1.781250 
O.OOOOOOOE+00 
0.7501E+08 
0.7330000E-03 
O.OOOOOOOE+00 
0 
0 
2 

0.3000000E+08 

0.3000000 

0.7330000E-03 

0.0 

1 

1 

N 

Y 

0 

0.3000000E+08 

0.3000000 

0.7330000E-03 

Y 

Y 

5.625000 

0.2500000 

1.250000 
N 

N 

H 

H 

10 

3 

1 

16.81250 

0.0 

0. OOOOOOOE+00 
0.0 



$ El * RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 70)( 2) 

$ E2 s AXIAL COMPONENT OF RING ECCENTRICITYC 2) 

$ GJ = TORSIONAL RIGIDITYC 2) 

$ RM * RING MATERIAL DENSITY (E.G. ALUMINUM* . 0002535) ( 2) 

$ E * YOUNG’S MODULUS OF RINGC 3) 

$ A * CROSS SECTION AREA OF RINGC 3) 

$ lY * MOMENT OF INERTIA ABOUT Y-AXIS (SEE FIG. ON P.70)( 3) 

$ IX * MOMENT OF INERTIA ABOUT X-AXISC 3) 

$ IXY * PRODUCT OF INERTIAC 3) 

$ El * RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 70)( 3) 

$ E2 * AXIAL COMPONENT OF RING ECCENTRICITYC 3) 

$ GJ * TORSIONAL RIGIDITYC 3) 

$ RM * RING MATERIAL DENSITY (E.G. ALUMINUM* . 0002535) ( 3) 

$ E * YOUNG’S MODULUS OF RINGC A) 

$ A * CROSS SECTION AREA OF RINGC 4) 

$ lY = MOMENT OF INERTIA ABOUT Y-AXIS (SEE FIG. ON P.70)( 4) 

$ IX * MOMENT OF INERTIA ABOUT X-AXISC 4) 

$ IXY * PRODUCT OF INERTIAC 4) 

$ El * RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 70) ( 4) 

$ E2 * AXIAL COMPONENT OF RING ECCENTRICITYC 4) 

$ GJ * TORSIONAL RIGIDITYC 4) 

$ RM = RING MATERIAL DENSITY (E.G. ALUMINUM* . 0002535) ( 4) 

$ E * YOUNG’S MODULUS OF RINGC 5) 

$ A * CROSS SECTION AREA OF RINGC 5) 

$ lY = MOMENT OF INERTIA ABOUT Y-AXIS (SEE FIG. ON P.70)( 5) 

$ IX * MOMENT OF INERTIA ABOUT X-AXISC 5) 

$ IXY = PRODUCT OF INERTIAC 5) 

$ El = RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 70)( 5) 

$ E2 * AXIAL COMPONENT OF RING ECCENTRICITYC 5) 

$ GJ = TORSIONAL RIGIDITYC 5) 

$ RM * RING MATERIAL DENSITY (E.G. ALUMINUM* . 00 02535 ) C 5) 

$ E * YOUNG’S MODULUS OF RINGC 6) 

$ A * CROSS SECTION AREA OF RINGC 6) 

$ lY = MOMENT OF INERTIA ABOUT Y-AXIS (SEE FIG. ON P.70)( 6) 

$ IX * MOMENT OF INERTIA ABOUT X-AXISC 6) 

$ IXY * PRODUCT OF INERTIAC 6) 

$ El * RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 70)( 6) 

$ E2 * AXIAL COMPONENT OF RING ECCENTRICITYC 6) 

$ GJ * TORSIONAL RIGIDITYC 6) 

$ RM * RING MATERIAL DENSITY (E.G. ALUMINUM* . 0002535) C 6) 

$ K*ELASTIC FOUNDATION MODULUS (E.G. LB/INX3{3)IN THIS SEG. 

$ LINTYP* INDICATOR CO, 1, 2 OR 3) FOR TYPE OF LINE LOADS 
$ NLTYPE*CONTROL (0, 1,2,3) FOR TYPE OF SURFACE LOADING 
$ NWALL*INDEX Cl, 2, 4, 5, 6, 7, 8) FOR WALL CONSTRUCTION 
$ E * YOUNG’S MODULUS FOR SKIN 

$ U = POISSON’S RATIO FOR SKIN 

$ SM *MASS DENSITY OF SKIN (E.G. ALUM. =.00025 LB-SECXX2/INXX4) 

$ ALPHA = COEFFICIENT OF THERMAL EXPANSION 

$ ANRS * CONTROL CO OR 1) FOR ADDITION OF SMEARED STIFFENERS 

$ SUR * CONTROL FOR THICKNESS INPUT (0 OR 1 OR -1) 

$ ARE THERE STRINGERS (PLEASE ANSWER Y OR N)? 

$ ARE THERE RINGS (PLEASE ANSWER Y OR N)? 

$ K2 =CONTROL (0 OR 1) FOR INTERNAL OR EXTERNAL RINGS 
$ E2 * RING MODULUS 

$ U2 * RING POISSON RATIO 

$ RGMD* RING MASS DENSITY 

$ IS THE RING CROSS SECTION CONSTANT IN THIS SEGMENT? 

$ IS THE RING CROSS SECTION RECTANGULAR (Y OR N)? 

$ D2 = ARC LENGTH BETWEEN ADJACENT RINGS (CONSTANT) 

$ T2 * THICKNESS OF RING (CONSTANT) 

$ H2 = HEIGHT OF RING (CONSTANT) 

$ DO YOU WANT TO PRINT OUT THE C(I,J) AT MERIDIONAL STATIONS? 

$ DO YOU WANT TO PRINT OUT DISTRIBUTED LOADS ALONG MERIDIAN? 

$ 

$ SEGMENT NUMBER 33333333 
$ NMESH * NUMBER OF NODE POINTS C5 = MIN.; 98 * MAX.)( 3) 

$ NTYPEH* CONTROL INTEGER (1 OR 2 OR 3) FOR NODAL POINT SPACING 

$ NSHAPE* INDICATOR (1,2 OR 4) FOR GEOMETRY OF MERIDIAN 
$ R1 = RADIUS AT BEGINNING OF SEGMENT (SEE P. 66) 

$ Z1 * AXIAL COORDINATE AT BEGINNING OF SEGMENT 

$ R2 * RADIUS AT END OF SEGMENT 
$ Z2 * AXIAL COORDINATE AT END OF SEGMENT 
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FILE: BOSSO DATA 



A1 



0 

3 

2.500000 

N 

0 

0.0 

0 

0 

2 

O.IOSOOOOE-HQS 

0.3200000 

0.2535000E-03 

0.0 

0 

1 

N 

N 

H 

H 

-1 

0 

0 

0 

1 

30 

O.OOOOOOOE+00 

O.OOOOOOOE+00 

O.OOOOOOOE+00 

1 



IMP = INDICATOR FOR IMPERFECTION (0=NONE» I=SOME) 

NTYPEZ= CONTROL (1 OR 3) FOR REFERENCE SURFACE LOCATION 
ZVAL = DISTANCE FROM LEFTMOST SURF. TO REFERENCE SURF. 

DO YOU WANT TO PRINT OUT R(S), RUS). ETC. FOR THIS SEGMENT? 
NRINGS= NUMBER (MAX=20) OF OISCRETE RINGS IN THIS SEGMENT 
K=ELASTIC FOUNOATION MOOULUS (E.G. LB/IN5tX3)IN THIS SEG. 
LINTYP= INOICATOR (0» 1, 2 OR 3) FOR TYPE OF LINE LOAOS 
NLTYPE=CONTROL (0,1»2,3) FOR TYPE OF SURFACE LOAOING 
NWALL=INDEX (1, 2» 4, 5, 6, 7, 3) FOR WALL CONSTRUCTION 
E = YOUNG'S MODULUS FOR SKIN 

U = POISSON'S RATIO FOR SKIN 

SM =MASS OENSITY OF SKIN (E.G. ALUM.*. 00025 LB-SEC3«x2/IN5tx4) 
ALPHA = COEFFICIENT OF THERMAL EXPANSION 

ANRS 5 CONTROL ( 0 OR 1) FOR AOOITION OF SMEARED STIFFENERS 
SUR = CONTROL FOR THICKNESS INPUT (0 OR 1 OR -1) 

00 YOU WANT TO PRINT OUT THE C(I,J) AT MERIDIONAL STATIONS? 
DO YOU WANT TO PRINT OUT DISTRIBUTED LOADS ALONG MERIDIAN? 

GLOBAL DATA 

NLAST = PLOT OPTIONS (-1=N0NE, 0=GEOMETRY, 1=U,V,W) 

NOB s STARTING NUMBER OF CIRC. WAVES (BUCKLING ANALYSIS) 
NMINB s MINIMUM NUMBER OF CIRC. WAVES (BUCKLING ANALYSIS) 

NMAXB = MAXIMUM NUMBER OF CIRC. WAVES (BUCKLING ANALYSIS) 

INCRB = INCREMENT IN NUMBER OF CIRC. WAVES (BUCKLING) 

NVEC 5 NUMBER OF EIGENVALUES FOR EACH WAVE NUMBER 
P = PRESSURE OR SURFACE TRACTION MULTIPLIER 

TEMP = TEMPERATURE RISE MULTIPLIER 

OMEGA = ANGULAR VEL . ABOUT AXIS OF REVOLUTION (RAD/SEC) 
NUMBER OF POLES (PLACES WHERE R=0) IN SEGMENT( 1) 

CONSTRAINT CONDITIONS FOR SEGMENT NO. I 1 1 1 

IPOLE = NOOAL POINT NUMBER OF POLE, IPOLE( 1) 

AT HOW MANY STATIONS IS THIS SEGMENT CONSTRAINED TO GROUND? 
IS THIS SEGMENT JOINED TO ANY LOWER-NUMBERED SEGMENTS? 
NUMBER OF POLES (PLACES WHERE R=0) IN SEGMENT( 2) 



1 
1 
1 

10 

1 

1 

1 

1 

O.OOOOOOOE't'OO 

O.OOOOOOOE+OO 

Y 

1 

H 

H 

10 

0 

Y 

1 

1 

2 
45 

1 

1 

1 

1 

O.OOOOOOOE+00 

O.OOOOOOOE+00 

Y 

Y 

Y 

Y 
N 



CONSTRAINT CONDITIONS FOR SEGMENT NO. 2 2 2 2 

AT HOW MANY STATIONS IS THIS SEGMENT CONSTRAINED TO GROUND? 

$ IS THIS SEGMENT JOINED TO ANY LOWER-NUMBERED SEGMENTS? 

$ AT HOW MAY STATIONS IS THIS SEGMENT JOINED TO PREVIOUS SEGS.? 
INODE = NODE IN CURRENT SEGMENT (ISEG) OF JUNCTION, INOOE( 1) 
JSEG = SEGMENT NO. OF PREVIOUS SEGMENT INVOLVED IN JUNCTION 
JNOOE = NOOE IN PREVIOUS SEGMNT (JSEG) OF JUNCTION 
IUSTAR= AXIAL OISPLACEMENT (0=NOT SLAVED, 1=SLAVE0) 

IVSTAR= CIRCUMFERENTIAL OISPLACEMENT (0=N0T SLAVED, 1=SLAVED) 
IWSTARs RAOIAL OISPLACEMENT (0=NOT SLAVED, I=SLAVED) 

ICHI = MERIDIONAL ROTATION (0=NOT SLAVED, 1=SLAVED) 

01 = RAOIAL COMPONENT OF JUNCTURE GAP 

02 = AXIAL COMPONENT OF JUNCTURE GAP 

IS THIS CONSTRAINT THE SAME FOR BOTH PREBUCKLING AND BUCKLING? 
NUMBER OF POLES (PLACES- WHERE R=0) IN SEGMENT( 3) 



CONSTRAINT CONDITIONS FOR SEGMENT NO. 3 3 3 3 

IPOLE = NODAL POINT NUMBER OF POLE, IPOLE( 1) 

AT HOW MANY STATIONS IS THIS SEGMENT CONSTRAINED TO GROUND? 

IS THIS SEGMENT JOINED TO ANY LOWER-NUMBERED SEGMENTS? 

AT HOW MAY STATIONS IS THIS SEGMENT JOINED TO PREVIOUS SEGS.? 
INOOE = NOOE IN CURRENT SEGMENT (ISEG) OF JUNCTION, INOOE( 1) 
JSEG = SEGMENT NO. OF PREVIOUS SEGMENT INVOLVED IN JUNCTION 
JNOOE = NOOE IN PREVIOUS SEGMNT (JSEG) OF JUNCTION 
IUSTAR= AXIAL DISPLACEMENT (0=NOT SLAVED, 1=SLAVED) 

IVSTAR= CIRCUMFERENTIAL OISPLACEMENT (0=NOT SLAVED, 1=SLAVED) 
IWSTARs RADIAL DISPLACEMENT (0=NOT SLAVED, 1=SLAVED) 

ICHI = MERIDIONAL ROTATION (0=NOT SLAVED, 1-SLAVED) 

D1 = RAOIAL COMPONENT OF JUNCTURE GAP 
02 = AXIAL COMPONENT OF JUNCTURE GAP 

IS THIS CONSTRAINT THE SAME FOR BOTH PREBUCKLING ANO BUCKLING? 
DO YOU WANT TO LIST OUTPUT FOR SEGMENT! 1) 

00 YOU WANT TO LIST OUTPUT FOR SEGMENT! 2) 

DO YOU WANT TO LIST OUTPUT FOR SEGMENT! 3) 

00 YOU WANT TO LIST PREBUCKLING RESULTANTS ANO RING FORCES? 
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2. ACESNID INPUT DATA 

The ACESNID computer code determines the virtual mass 
array, a quantity producing the late-time contribution of 
the DAA. Only one ACESNID calculation is required for all 
circumferential harmonics of the shell. The bands dividing 
the areas of the shell are identified in the NLEFT , NCYL, 
and NRITE entries, while the fluid properties are specified 
in the RHOFL and VSOUND row. 



FILE-. ACEl DATA A1 



VIRTUAL MASS FOR RING-STIFFENED CYLINDER WITH FLAT ENDS, N=0,1 

01110/ NSTART,NFINIS,NFREQ,NVMASS,NCHECK 

15 51 15 3 0 / NLEFT, NCYL, NRITE, NWBOSG.NSYMF 

1 6 10 0 0 / NORDER,NFENDS,NFCENT,NFCMPT,NOMIT 

9.59684E-5 5.833E4 0.005 / RHOFL , VSOUND, ERR 

1111110111/ OUTPUT FLAGS 

1.0 / CPS(l) 

1 1 1 12 / SEFS, JSGBEG, JPTBEG, JSGEND, JPTEND,LEFT 

2 1 2 47 / SEFS, JSGBEG, JPTBEG, JSGEND, JPTEND, CENTRAL 

3 1 3 12 / SEFS, JSGBEG, JPTBEG, JSGEND, JPTEND, RIGHT 

1 1 1 12 / SOURCES, JSGBEG, JPTBEG, JSGEND, JPTEND, LEFT 

2 1 2 47 / SOURCES, JSGBEG, JPTBEG, JSGEND, JPTEND, CENTRAL 
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3. PIFLASH INPUT DATA 

The PIFLASH code takes data from the shell mode files 
from B0S0R4 and the virtual mass file from ACESNID and 
combines them to create a "shell-fluid file" containing all 
the information concerning the shell and its acoustic medium 
for further calculations. Since only the breathing modes 
are activated in end-on loading, all torsional modes are 
dropped from the shell file. This is accomplished in the 
NJUSE rows where the modes of Table I are identified. 



FILE: PIEl DATA A1 



20000/ NUMBER, NTORSN,NPTM,NSYMS,NSYMP 
3333/ (NWETSG(K),NUSESG(K),K=1,2) 

386.4 / GRAVITY 

SHELLOO / SMF, N=0 

SHELLOl / SMF, N=1 

19 2 / (NJUSE(J), J=1,NITEMS) 

123/ (KORSG(K),K=1,NKORSG) 

2 4 6 7 9 12 13 14 19 20 21 22 23 25 26 27 28 29 30 / JUSE, N=0 
12/ JUSE, N=1 
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4. 



SAPIV INPUT DATA 



The substructure is modeled using the standard SAPIV 
code with slight modifications by the developers of ELSHOK. 
The diaphragm is modeled using twenty-four plates and 
twenty-five nodes. The location of each node is identified 
in cartesian coordinates. Each plate is 1.06 in thick for 
the resonant case. Nodes eighteen through twenty-five are 
rigidly attached to the shell in the last eight lines of the 
code. The SAPIV code determines the fixed-base modes and 
corresponding natural frequencies as well as the uncon- 
strained mass and stiffness matrices for the substructure. 
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FILE: SAP23 



DATA 



A1 



IDEALIZED INTERNAL EQUIPMENT AT FRAME 9 
25 10 10 100000/ NUMNP,NELTYP, . . . ,N10SV,NRIGID 
0 1 0 1 1 1 1 1 0.0 0.0 0.0 0 0.0 
0 2 0 0 0 1 0 0 0.0 0.0 5.60A 0 0.0 
03000100 0.0 3.963 3.963 0 0.0 
OAOOOIOOO.O 5.60A 0.0 0 0.0 
05000100 0.0 3.963 -3.963 0 0.0 
0 6 0 0 0 1 0 0 0.0 0.0 -5.604 0 0.0 
07000100 0.0 -3.963 -3.963 0 0.0 
08000100 0.0 -5.604 0.0 0 0.0 
09000100 0.0 -3.963 3.963 0 0.0 



0 


10 


0 


0 


Q 


1 


0 


0 


0.0 


0.0 11.208 0 0.0 


0 


11 


0 


0 


0 


1 


0 


0 


0.0 


7.925 7.925 0 0.0 


0 


12 


0 


0 


Q 


1 


0 


0 


0.0 


11.208 0.0 0 0.0 


0 


13 


0 


0 


0 


1 


0 


0 


0.0 


7.925 -7.925 0 0.0 


0 


14 


0 


0 


0 


1 


0 


0 


0.0 


0.0 -11.208 0 0.0 


0 


15 


0 


0 


0 


1 


0 


0 


0.0 


-7.925 -7.925 0 0.0 


0 


16 


0 


0 


0 


1 


0 


0 


0.0 


-11.208 0.0 0 0.0 


0 


17 


0 


0 


0 


1 


0 


0 


0.0 


-7.925 7.925 0 0.0 


0 


18 


0 


0 


0 


1 


Q 


0 


0.0 


0.0 16.813 0 0.0 


0 


19 


0 


0 


0 


1 


0 


0 


0.0 


11.888 11.888 0 0.0 


0 


20 


0 


0 


0 


1 


0 


0 


0.0 


16.813 0.0 0 0.0 


0 


21 


0 


0 


0 


1 


0 


0 


0.0 


11.888 -11.888 0 0.0 


0 


22 


0 


0 


Q 


1 


0 


0 


0.0 


0.0 -16.813 0 0.0 


0 


23 


0 


0 


0 


1 


0 


0 


0.0 


-11.888 -11.888 0 0.0 


0 


24 


0 


0 


0 


1 


0 


0 


0.0 


-16.813 0.0 0 0.0 


0 


25 


0 


0 


0 


1 


0 


0 


0.0 


-11.888 11.888 0 0.0 



6 24 1 / TYPE 6 (PLATE AND SHELL ELEMENTS) 

1 .733E-03 0.0 0.0 0.0 

'32.967E+06 9.89E+06 0.0 32.967E+06 0.0 11.538E+06 

0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 / ELEMENT LOAD MULTIPLIERS 

118900101. 06 000 

219200101. 06 000 
31230010 1. 06 000 

413400101. 06 000 

514500101. 06 000 
615600 101. 06 000 
71670010 1.06 000 

817800101. 06 000 

9 2 10 11 3 0 1 1 1.06 000 

15 8 16 17 9 0 1 0 1.06 0 0 0 

16 9 17 10 2 0 1 0 1.06 000 

17 10 18 19 11 0 1 1 1.06 0 0 0 

23 16 24 25 17 0 1 0 1.06 0 0 0 

24 17 25 18 10 0 1 0 1.06 0 0 0 

0 0 0.0 0.0 0.0 0.0 0.0 0.0 / CONC MASSES 

0.0 0.0 0.0 0.0 / ELEMENT LOAD MULTIPLIERS 
0 0 0 0.0 0.0 0 0.0 / IFPR,IFSS, . . .,NFO,SHIFT 
80/ NFACE,NSYMTY 

18 1 1 1 1 0 0 

19 1 1 1 1 0 0 

20 1 1 1 1 0 0 

21 1 1 1 1 0 0 

22 1 1 1 1 0 0 

23 1 1 1 1 0 0 

24 1 1 1 1 0 0 

25 1 1 1 1 0 0 

END OF DATA (DUMMY TITLE) 
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5 



PICRUST INPUT DATA 



The PICRUST code takes data from the substructure mode 
file from SAPIV and reorganizes it to facilitate the 
solution of the transient response equations. The location 
of the substructure attachment points to the shell is 
identified for each connecting node. 



FILE: PICRUSTl DATA A1 



1111011100 

0111110101/ OUTPUT FLAGS 1-20 
200/ NJUSE,NHWS0B,NHWBAR 
18 2 11 0.0 

19 2 11 45.0 

20 2 11 90.0 

21 2 11 135.0 

22 2 11 180.0 

23 2 11 225.0 

24 2 11 270.0 

25 2 11 315.0 / NIPSUB,LB0SEG,LB0SPT,ANGDEG 
0.0 / DEGROT 
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6. USLOB INPUT DATA 

The two USLOB input codes for the taper and conventional 
charges appear on the following page. USL46 DATA is the 
code for the taper charge. It contains nine discrete data 
points from the pressure-t ime history in seconds and pounds 
per square inch. There are 1600 time points used in the 
numerical integration of the governing equations with 
information saved every sixteenth point, so only 100 
velocity-time data points are retained for EASYPLOT display. 
The charge is placed along the negative X-axis at 840 
inches. USL47 DATA contains the information for the 
conventional charge. The only difference between the two 
codes is a charge identification entry (NCHRG) and the 
charge weight and constants specified for the conventional 
charge instead of the pressure-t ime history (WCHRG-THEXP) . 
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FILE: USL«»6 DATA A1 



1599 16 2 9 2 2 1 / NTIME, NSKIP , NCHRG, NQUAD, NFINE, KOUPLE, NSUBS 
5.0E-05 -340.0 0.0 0.0 / DELT,XLOAD, RLOAD,SURCUT 
9 / NSHAPE 

0.0 240.0 0.0002 320.0 0.00053 350.0 0.00243 305.0 0.00294 160.0 

0.0034 115.0 0.0044 70.0 0.0054 48.0 0.00565 0.0 / (TC( J ), PCCJ ), J=1 , NSHAPE) 

0.0 / DECAY 

0100000000 

0010000000/ OUTPUT FLAGS 1-20 
4 / NPTSHL 

1 1 2 11 2 46 3 12 / LBOSEG,LBOSPT 
1 / NPTSUB 



FILE: USL47 DATA A1 



1599 16 1 9 2 2 1 / NTIME, NSKIP, NCHRG, NQUAD, NFINE, KOUPLE, NSUBS 
5.0E-05 -840.0 0.0 0.0 / DELT , XLOAD, RLOAD, SURCUT 

352.0 3.3354314E+5 1.144 3.03131E-5 -.247 / WCHRG, PZMLT,PZEXP,THMLT,THEXP 
0100000000 

0010000000/ OUTPUT FLAGS 1-20 
4 / NPTSHL 

1 1 2 11 2 46 3 12 / LBOSEG,LBOSPT 
1 / NPTSUB 
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7. PUSLOB INPUT DATA 

The PUSLOB code produces TEKTRONIX plots of the 
velocity-time histories for specified points on the shell or 
substructure. It also generates velocity punch card files 
which are converted to a format for EASYPLOT using the 
PLOTCONVl code. 
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FILE: PUSL^6 DATA 



A1 



1599 16 1 1 1 / NTIME,NSKIP,NSUBS,NTEK,NCARD 
5.0E-05 1000.0 1.0 / DELT,XMULT,YMULT 
VELOCITY PROFILE 
TIME (MSEC) 

VELOCITY (IN/SEC) 

NPTSHL-W 
-FWD 1/01 
1 1 -1 1 0.0 
V-FWD 3/12 

3 12 -1 1 0.0 / LBOSEG, LBOSPT,NSPHC,NANG,ANGDEG 
13 NPTSHL-U 
V-FWD 1/01 
1 1 1 1 0.0 • 

V-FWD 2/02 
2 2 0 1 0.0 
V-FWD 2/07 
2 7 0 1 0.0 
V-FWD 2/11 
2 11 0 1 0.0 
V-FWD 2/15 
• 2 15 0 1 0.0 

V-FWD 2/20 
2 20 0 1 0.0 
V-FWD 2/2A 
2 2A 0 1 0.0 
V-FWD 2/28 
2 28 0 1 0.0 
V-FWD 2/33 
2 33 0 1 0.0 
V-FWD 2/37 
2 37 0 1 0.0 
V-FWD 2/Al 

2 <4l 0 1 0.0 

V-FWD 2/A6 

2 46 0 1 0.0 
V-FWD 3/12 

3 12 1 1 0.0 / LBOSEG, LBOSPT,NSPHC, NANG, ANGDEG 

2 NPTSHL-V 

V-FWD 1/01 
11-11 0.0 
V-FWD 3/12 

3 12 -1 1 0.0 / LBOSEG, LBOSPT,NSPHC, NANG, ANGDEG 
1 NPTSUB 

PLT CNTR 
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